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Abstract. - We study the effect of weak and dilute disorder on the order parameter equation and 
transition temperature of a Pomeranchuk-type Fermi-surface instability using replica mean field 
theory. We consider the example of a phase transition to a d x 2_ y 2 type Fermi surface distortion, 
and show that, in the regime where such a transition is second order, the transition temperature 
is reduced by disorder in essentially the same way as that for a d-wave superconductor. We argue 
that observing this disorder dependence of metal-to-metal transition is a useful indicator of a finite 
angular momentum Fermi surface distortion. 



I 

C 

O 
o 



> 

in 
in 

ON 

^6 
o 
i> 
o 



Introduction. — In recent years a growing number 
of correlated systems have been found to exhibit thermo- 
dynamic phase transitions between metallic states. Ex- 
amples include the 17K transition in URu2Si2 [1] and 
the transitions around the metamagnetic quantum criti- 
cal endpoint of Sr 3 Ru207 [2]. In contrast to the more fa- 
miliar superconducting or magnetic instabilities, the order 
parameter which presumably develops at these transitions 
appears to be transparent or only weakly coupled to most 
experimental probes. Thus identifying the nature of this 
"dark" or "hidden" order is a challenging problem. Many 
years ago Pomeranchuk [3] found a condition for instabili- 
ties between metallic states characterized by Fermi surface 
shape distortions. It has been argued that this instability 
occurs in quantum Hall systems [4, 5] and is the origin of 
the transitions in URu2Si2 [6] and Sr 3 Ru207 [2]. How- 
ever, the key question remains how to identify this sort of 
order-particularly if bulk changes are masked by domains 
formation. In this Letter we calculate the form of the dis- 
order dependence of the transition temperature and find 
it to have a characteristic signature of momentum space 
distortions of the metallic Fermi surface. 

We are motivated by the empirical similarity in the 
way the mysterious phase in the bilayer ruthenate, 
Sr3Ru2C>7 [2], and the superconductivity in the related 
single-layer compound, Sr2Ru04 [7], are both eliminated 
with very low levels of disorder. In the latter case quan- 
titative comparison of the strong disorder dependence of 
the superconducting transition temperature T c to the well- 
known form [8] has become the de-facto signature of non- 



zero orbital angular momentum pairing (in the absence of 
phase sensitive methods). In contrast s-wave supercon- 
ductors are insensitive to disorder [9]. Could a similar 
dependence be used to diagnose the Pomeranchuk transi- 
tion? 

In this Letter we show how the Pomeranchuk instabil- 
ity is an analogue of non s-wave superconductivity but in 
the particle-hole rather than particle-particle channel. We 
exploit this using a combination of standard methods to 
show that, not only might one similarly expect a sensitiv- 
ity to disorder but that the precise form of the disorder de- 
pendence of T c is, under certain circumstances, identical to 
that of unconventional superconductors. This provides a 
quantitative test of the "dark order" metallic phase which 
parallels that now used for unconventional superconduc- 
tivity. Moreover, since we show that the Pomeranchuk 
instability is strongly suppressed in presence of weak dis- 
order, our results may provide an explanation as to why 
this rather subtle metal to metal transition is not observed 
more often in nature. 

Model. — To study effects of disorder in a simple 
model of a Pomeranchuk instability, we consider elec- 
trons on a two dimensional (2D) tight-binding lattice with 
a quadrupolar interaction that has been studied exten- 
sively [10-13] 
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/ k — [-0 k | "0kjj i s tne spinor creation operator. 
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dcncc: Vq(k, p) = <7 q k p , where the d x 2_ y 2 form factor 
<Pk = cosk x — cosfcj,. Kee et al. [11] found that the d xy 
component of the quadrupolar interaction usually does not 
acquire an expectation value, which we thus drop. 

We add a weak, dilute disorder potential that couples to 
the electron density i?dis = / d d x£(x) tp^ (x)i(>(x) assum- 
ing, for simplicity, delta-correlated, static (quenched) dis- 
order Pr(£) oc exp-J d d x S > 2 (x)/2D. Here D = \/2itN t 
is the strength of the disorder potential, with N the den- 
sity of state at the Fermi surface, and t is the disorder 
scattering time. We do not here consider stronger or cor- 
related disorder, where the disorder potential may couple 
directly to the Fermi surface distortion (Pomcranchuk) or- 
der parameter [14]. 

Methods. — For quenched disorder, one needs to 
disorder-average the free energy instead of the partition 
function. One standard method that works also for in- 
teracting systems is the replica trick [15,16], based on 
the identity: InZ — lim„^o(^" — l)/ n • Note that we 
have also derived the results presented below using dia- 
grammatic perturbation theory [17] but we found that the 
replica method makes the parallel with unconventional su- 
perconductivity explicit. The idea is to replicate n copies 
of the partition function Z, disorder-average Z n , and fi- 
nally take the limit n — > to get the disorder-averaged 
free energy. Since we have taken a simple Gaussian dis- 
tribution for the disorder potential, this disorder-average 
is readily done to give the disorder-induced interaction 4- 
fermion term (in momentum representation) 
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where the subscript n, m on the electron operators denote 
the Matsubara frequencies w n ,w m , also J2 n m refers to 
Matsubara frequencies summation, and a,(3 — 1, . . . n are 
replica indices. 

To derive a low energy effective theory we follow Be- 
litz and Kirkpatrick [16] and consider the disorder-induced 
interaction with all momenta near the Fermi surface. 
There are three possible ways of pairing up the scatter- 
ing: (1) is the small angle (or direct) scattering with 
k 2 = ki + q, (2) is the large angle (or exchange) scat- 
tering with k 4 — ki + q, and (3) is the pair (or 2kp) 
scattering where k3 = ki + q. The momentum transfer q 
is now restricted to be small (with a cut-off much smaller 
than the Fermi momentum). It can be shown that type 
(1) only leads to a renormalization of the chemical poten- 
tial and we drop it from now on. Type (3) couples to the 
superconducting order parameter, but not to the Pomer- 
anchuk one, and we can show that this term does not have 
any effect on the Pomeranchuk order at mean field level 
so we neglect it. However, for a superconductor, the type 
(3) term generates a vertex correction that for an s-wave 
superconductor cancels the propagator correction due to 



type (2) term, thereby rendering it insensitive to disorder 
(Anderson's theorem) [17]. Thus, the disorder-induced in- 
teraction is 
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The full low energy effective action after disorder- 



averaging is thus S = Yl a (^0 + ^ir 
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To decouple the four fcrmion interaction terms, we 
introduce the Q matrix via essentially a Hubbard- 
Stratonovich transformation (generalizing Ref. [18]) 
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where i,j label the spinor components. 

Since the replica-Q-matrix method has already been 
comprehensively reviewed in Ref. [16], [18], we here only 
sketch out its application to the Pomeranchuk instabil- 
ity in the presence of weak quenched disorder. Assuming 
that at the saddle point, there is replica symmetry and 
spin symmetry, the homogeneous and un-retarded ansatz 
for the saddle-point of this action is 
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Note that only one function Q„k is needed here (unlike 
for magnets or superconductors), because both disorder 
and the quadrupolar interaction induces a self-energy iA„k 
that enters in the same way into the propagator renormal- 
ization. A„k is the Fourier transform dual of Q n v. [18] and 
thus has the same structure as Q„k- 

With the ansatz eq. 7, the saddle point action becomes 
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Note that only the q = component of Vq(k, p), i.e. 
Vo(k, p) matters, because of the assumption of spatial ho- 
mogeneity in the saddle point ansatz Q nk . Also we have 
assumed replica symmetry: we can a posteriori justify 
this by noting that as we shall see, there are no indi- 
cations in the free energy of further instabilities in the 
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Fermi surface distorted phase, unlike in the classical spin 
glass case which does demand replica symmetry breaking. 
Presumably this has to do with the much simpler free en- 
ergy landscape in the Pomeranchuk case, indicating lack of 
glassiness in our system. Replica symmetry means we can 
drop the replica indices from now on. The saddle point 
equations SSgp/SQnk = and (55 sp /<5A nk = give 
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First, lets check that the ansatz eq. 7 recovers known 
results. For free electrons with quenched disorder, setting 
Vo(k, p) = leads to the standard Born approximation 
result; at the saddle point, Q„ k is just the electron prop- 
agator with a disorder-induced lifetime r: Q nk ~ Gnk = 
[iu> n — £k + ^sgn^n)] -1 - For clean electrons with the 
quadrupolar interaction, setting r — > oo, and assuming 
a spatially homogeneous order parameter (i.e., only the 
q = component of V^(k, p) is involved), we define the 
Pomeranchuk order parameter 

^ u A =T^Fo(k,p)(^ p ^p) = 20 k T^#pQ„ p . 
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We then recover the clean case mean field order parameter 
equation A = 2g J^u^ h («* + </>kA ), where f T (x) = 
[expa;/r+ is the usual Fermi distribution. 

Results. — Now we consider the case of weakly dis- 
ordered electrons with an interaction favoring a Pomer- 
anchuk instability, eqs. 11,9,10 lead to the order parame- 
ter equation 
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This is in fact the non s-wave, non-magnetic analogue of 
the Stoner instability of a ferromagnet. The extra angular 
dependence in the momentum sum means that this order 
parameter equation does not reduce to the clean case. (By 
contrast, the s-wave Stoner Pomeranchuk instability is un- 
affected, to leading order, by impurities [19]). Thus, just 
as for non s-wave superconductors, there is no Anderson's 
theorem for / ^ Pomeranchuk instabilities, because of its 
angular dependence in momentum space, eq. 12 simplifies 
to: 
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As expected, disorder smears out the Fermi distribution 
to give the digamma function ip(x): crudely speaking, dis- 
order raises the effective temperature. 

In contrast to weak-coupling superconductivity where 
there is always a second-order transition for arbitrary 
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Fig. 1: The effect of disorder on the magnitude of the order pa- 
rameter A (all energies and scattering rates expressed in units 
of 2t). (a) At T = as disorder decreases the lifetime, r, the 
critical coupling, g/2ir 2 increases and a first order transition 
becomes second order, (b) At fixed coupling, the order param- 
eter and T c are rapidly suppressed by weak disorder. 



weak interactions, the Pomeranchuk instability requires 
a critical coupling and the transition can either be first 
or second-order. To check the order of the transition, we 
need to evaluate the free energy. Substituting the mean 
field equations 9,10 into the saddle point action eq. 8, to- 
gether with the approximation (just as for disordered free 
electrons case) ^ p Q n p ~ — wr./VoSgna; n and the definition 
of the order parameter parameter (eq. 11) , the saddle 
point free energy becomes 



F sp = -Trln 



—i u n + 
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where the renormalized dispersion e k = e k + /i + Ao</> k . 
Defining the renormalized density of state Na(c) = 
2 k S(e — e k ), the order parameter equation 12 becomes 



^4T/ 4t d e^Relnr 
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For quantitative results for effects of disorder, we eval- 
uate the free energy and mean field equations numeri- 
cally for the 2D square lattice with bare dispersion e k = 
— 2t (cos k x + cos ky) — /i, and the d-wave form factor: 
</> k = cos fez — cos k y , with the corresponding renormalized 
density of states (see Ref. [11, 12] for the actual form). 
In the following, energies and scattering rates, r _1 , are 
measured in units of 2t. 

First we consider how disorder changes the evolution of 
the order parameter. In Fig. 1(a) we see for fixed chemi- 
cal potential, /i, how the Pomeranchuk order parameter is 
modified at T = by scattering. Disorder both increases 
the critical coupling and turns at T — from first order in 
the clean limit, to second order. In Fig. 1(b) we see that 
rather weak disorder dramatically reduces T c . 

We next consider the disorder dependence of the Pomer- 
anchuk transition temperature, T c . In Fig. 2(a) we choose 
a parameter region where the clean transition is second 
order {g/2t2ir 2 = 0.051 and (i/2t = 0,0.05). Since the or- 
der parameter goes smoothly to zero, we can simplify the 
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Fig. 2: The Pomeranchuk transition temperature T c relative 
to the clean value T®, as a function of dimensionless disorder 
1/tT® at a parameter regime where the Pomeranchuk transi- 
tion is second order, (a) At y, — the numerical solution of the 
order parameter equation 14 (small points) coincides with the 
Abrikosov-Gorkov result of eq. 16 (solid line). The large data 
points show a fit to published data [2,21] on SraRu207 at 7.95T 
where the resistivity minimum and the zero-field residual resis- 
tivity indicated T c and 1/r respectively, (b) At fj, — 0.05 there 
are small deviations from the Abrikosov-Gorkov form. These 
are corrections from the weak-coupling approximation and are 
related to the detailed structure of the density of states we 
used. 

order parameter equation (eq. 12) to get the second order 
transition order parameter equation 

-- = TV^- =— -. (15) 

3 f^f k K+27Sgn^)2 + 4 

We then note that this order parameter equation is iden- 
tical to the one determining the critical temperature for 
(i-wave superconductor with non-magnetic disorder. Note 
that this is true only for the T c equation for a second or- 
der transition: there are first order transitions at larger /i, 
and furthermore, the full order parameter equation for the 
disordered Pomeranchuk instability has a different form to 
the d-wave superconductor with disorder. 

Thus, from eq. 15, we get the familiar Abrikosov- 
Gor'kov form [20] for the disordered gap equation 

In figs. 2 (a) and (b), the solid curve is this universal 
Abrikosov-Gorkov form, while the small data points are 
direct numerical evaluation of the general (i.e. not just for 
second order transition) order parameter equation 14. At 
n/2t = 0, the direct evaluation coincide with the universal 
form, while for finite chemical potential fx/2t — 0.05 shown 
in Fig. 2 (b), some small deviation can be seen at larger 
disorder, due to the approximation in the radial k-integral 
that goes into deriving eq. 16. 

We have also taken existing data of on Sr3Ru2 07 at 
7.95T [2,21] where the resistivity minimum and the zero- 
field residual resistivity are taken to indicate T c and 1/r 
respectively, and plotted them as large dots in Fig. 2(a). 
The reasonable fit shows that the putative transition in 




Fig. 3: Suppression of T c as a function of scattering rate in 
units of 2t for a Fermi surface away from the van Hove point. 
The transition is driven first order at low temperatures. 

Sr 3 Ru2C>7 (which experimentally, is found not to be of su- 
perconductivity type) does follow the universal Abrikosov- 
Gorkov form for disordered Pomeranchuk transition, even 
if our actual model interaction of eq. 1 may be too sim- 
plistic. 

Finally we consider the disorder dependence on the tran- 
sition temperature where, at low temperatures, the tran- 
sition can become first order [22] such as when the sys- 
tem is further away from half-filling, e.g. with ji — 0.10 
(Fig. 3). Then, with larger disorder (larger 1/r), the tran- 
sition turns from second to first order (bold line). Sur- 
prisingly, the effect of increasing disorder is opposite to 
increasing T. Higher T smears out the Fermi function and 
leads to a smaller order parameter and eventually a second 
order transition results. [12] Of course, it is possible that 
including Gaussian fluctuations around the saddle point 
may turn first order transitions into second order ones in 
the presence of disorder. Future work is needed to resolve 
this issue. What is clear is that even at the mean field 
level, there is a strong suppression of T c for the Pomer- 
anchuk instability with increasing disorder. 

In summary, we have calculated the strong dependence 
of the d-wave Fermi surface distortion transition on di- 
lute disorder and shown that the effect is reminiscent both 
qualitatively and quantitatively of the strong dependence 
on impurities of the order parameter and T c in non s-wave 
superconductors. Our results suggest that detailed dis- 
order dependence of "hidden order" transitions could be 
used to indicate Pomeranchuk type order just as is the case 
for low T c unconventional superconductors. The extreme 
sensitivity to disorder [2] of the low-temperature metal- 
to- metal transition in Sr 3 Ru2 07 is suggestive of this [see 
Fig. 2(a)]. Our theoretical results are intended to moti- 
vate a more comprehensive systematic experimental study 
to compare with our quantitative predictions for the pu- 
tative ordered state in Sr 3 Ru207. 
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